clear all
global  dtSpikes ANoutput ANprobRateOutput
BF= 5000;
MAPparamsName='Normal';
AN_spikesOrProbability= 'spikes';
probeFrequency= BF;
probedB= 10;
sampleRate= 100000; % Hz
dt= 1/sampleRate; % seconds
PSTHbinWidth= 0.0005;
probeDuration= .200;  % seconds
rampDuration= 0.001;
beginSilenceDuration= 0.01; %(s)
beginSilence= zeros(1, round(beginSilenceDuration* sampleRate));
endSilenceDuration= 0.01 ; %(s)
endSilence=   zeros(1, round(endSilenceDuration* sampleRate));
allFiberTypes= {'HSR'};
%HSR= 230e-6;
%MSR= 125e-6;
%LSR= 70e-6;
%paramChanges={['IHCpreSynapseParams.tauCa=230e-6;']};
paramChanges={['']};

% generate probe tone
time1= dt: dt: probeDuration;
amp= 10^(probedB/20)*28e-6;
tone= amp*sin(2*pi*probeFrequency*time1);
rampTime= dt:dt:rampDuration;
ramp= [0.5*(1+cos(2*pi*rampTime/(2*rampDuration)+pi)) ones(1,length(time1)-length(rampTime))];
tone= tone.*ramp;
tone= tone.*fliplr(ramp);
        
inputSignal= [beginSilence tone endSilence];
overallDuration= dt*length(inputSignal);

